Effect of disorder on the NMR relaxation rate in two-band superconductors 
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We calculate the effect of nonmagnetic impurity scattering on the spin-lattice relaxation rate in 
two-band superconductors with the s-wave pairing symmetry. It is found that for the interaction 
parameters appropriate for MgB2 the Hebel-Slichter peak is suppressed by disorder in the limit 
of small interband impurity scattering rate. In the limit of strong impurity scattering, when the 
gap functions in the two bands become nearly equal, the single-band behavior is recovered with a 
well-defined coherence peak just below the transition temperature. 

PACS numbers: 74.25.Nf, 74.70.Ad, 74.20.-z 
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I. INTRODUCTION 



The discovery of superconductivity in MgB2 Llll has 



lit 

revived the interest in two-band superconductors yj . 
Most experiments on this compound Q are consistent 
with two distinct superconducting gaps Ajr and Ao- , with 
Act/Ajt ~ 2.63 5]. The only exception seems to be the 
nuclear magnetic resonance (NMR) relaxation rate T^^ 
[g,, in the superconducting state. The NMR relaxation 
rate measurements on ^^B nucleus by Kotegawa et al. 
[^Q found a small Hebel-SHchter coherence peak Q just 
below the superconducting transition temperature T^, al- 
though there are some questions as to whether the peak 
is present at all 0,^3- Moreover, by assuming that the 
contributions to the relaxation rate from the two bands 
are simply additive Kotegawa et al. found 7] a better 
fit to their data on ^^B nuclei by assuming that only one 
band participates in the nuclear relaxation process. 

In a recent work we have examined the NMR relax- 
ation rate in both conventional and unconventional clean 
two-band superconductors. The essence of our findings 
is that if the Fermi contact interaction between the nu- 
cleus and the conduction electrons dominates the relax- 
ation mechanism, then the measurements of T^^ probe 
the properties of the electron subsystem which are lo- 
cal in real space and hence extremely non-local in the 
momentum space fQ. We showed that for a two-band 
superconductor this gives rise to inter-band interference 
terms in T^^ , in addition to direct contributions from 
each band. 

Here we consider the effect of nonmagnetic impurity 
scattering on T-T ^ for a two-band superconductor. It is 
well known [la ll4ll5lll6| that the interband scattering 
by nonmagnetic impurities reduces the transition tem- 
perature of an s-wave two-band superconductor. In the 
case of MgB2 it has been argued (131 that one requires de- 
fects which produce lattice distortions in order to achieve 
large enough interband scattering which would lead to a 
significant suppression of Tc. Experimentally, such reduc- 
tion of Tc was observed on samples of MgB2 irradiated by 
fast neutrons [TsL IT^ I20l l2ll|. However, there is a finite 
interband scattering by impurities even in unirradiated 
and undoped MgB2 as evidenced by the break-junction 
tunneling experiments whose interpretation was es- 



sentially based on the work of Schopohl and Scharnberg 
1^35 on the tunneling density of states of a disordered 
two-band superconductor. The interband scattering rate 
from the low gap band to the high gap band used to fit 
the tunneling data in was comparable to the tran- 
sition temperature of 39 K. Thus it seems reasonable to 
investigate the effect of disorder on T^^ over a wide range 
of scattering rates. 

The NMR relaxation rate in a superconductor is de- 
termined by both normal and anomalous local density of 
states. In the case of a two-band superconductor these 
quantities depend on the interband scattering 0,13 H^l 
and on temperature in a rather complex way and one 
cannot a priori guess the effect of disorder on the temper- 
ature dependence of Tj~^ except in the limits of a small 
and very large interband impurity scattering rate. In the 
former case the peaks in the partial densities of states 
are somewhat reduced and broadened by the interband 
scattering compared to the clean system and one would 
expect the Hebel-Slichter coherence peak in Tj~^ below 
to be reduced and broadened by disorder. In the limit of 
a very large interband scattering such that the gap func- 
tions in the two bands become nearly equal (the Ander- 
son limit [12), the difference between the local densities 
of states in the two bands has no consequences and one 
expects to find the temperature dependence of T^^ char- 
acteristic of a single-band superconductor, with the size 
of the Hebel-Slichter peak determined by the strength of 
the electron-phonon coupling in the system [23 . |26| . In 
the case of MgB2, which is a medium coupling super- 
conductor 1^, the Hebel-Slichter peak is expected to be 
quite large in the Anderson limit. Indeed we find that 
our detailed numerical calculations of T^^ over a wide 
range of the interband impurity scattering rates confirm 
such a trend. 

The rest of the article is organized as follows. In Sec. 
im we summarize the strong coupling theory of the NMR 
relaxation rate in a disordered s-wave two-band super- 
conductor assuming that the Fermi contact interaction 
between the nuclear spin and the conduction electrons 
provides the dominant relaxation mechanism. In Sec. 
mil we give the results of our numerical calculations us- 
ing the interaction parameters for MgB2 0, 0, [l^ and 
in Sec. IIVI we give conclusions. 
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II. THEORY 



The relaxation rate of a nuclear spin due to the hyper- 
fine contact interaction with the band electrons is given 

by na 



1 



R^^--'^ lim ^-^^(-"). 



(1) 



where J is the hyperfine coupling constant, luq is the 
NMR frequency, and K^_{u;o) is the analytic continu- 



ation of the Fourier transform Kiy^ 
time correlator 



of the imaginary 



K{t) = -((Tr{S+iO,-iT)S-iO,0)))), 



(2) 



averaged over the impurity configurations. Here 
S±{r, —it) = e^"'^ S±{r)e~^'''^ , is the electron Hamil- 
tonian, and 

S+{r) = i;\{r)4,^{r), S^{r) = 4,\{r)4,^{r) (3) 



only the ladder impurity diagrams since the self-energies 
are calculated in the self-consistent second Born approx- 
imation. We have shown recently that in the case of 
a single three-dimensional band the contribution to the 
vertex function T from the ladder impurity diagrams is 
of the order 7/e_F, where 7 is the impurity scattering rate 
and Ei? is the Fermi energy, and that the impurity vertex 
corrections can therefore be neglected. The same applies 
to the contribution to from the ladder impurity dia- 
grams since the structure of the impurity ladder diagrams 
for the Fy S is completely analogous to the single-band 
case. Hence, we replace Fy in Eq. |(SJ) with matrix iT2 in 
computing the ratio of the spin-lattice relaxation rates in 
the superconducting and normal states. 

Next, by introducing the spectral representation for 
Gi{k,LUn) one can calculate the Matsubara sums in Eq. 

and then analytically continue the result to just above 
the real frequency axis ii^m (^0 + fl 1] . In the limit 
ujQ we obtain 



with tpl.{r) and tp^^{r) being the electron field operators 
{h — kB — 1 in our units, and we consider a system 
of unit volume) . contains both electron- phonon and 
screened Coulomb interactions as well as the scattering 
off randomly located nonmagnetic impurities. 

We assume that there are two spin-degenerate elec- 
tron bands in the crystal (the generalization to an arbi- 
trary number of bands is straightforward), and neglect 
the spin-orbit coupling. Then the spin operators ||2J) can 
be written in the band representation, using 



Rn 



+00 



where the densities of states (both normal and anoma- 
lous) are defined by 



(4) 



where Ui^kif) are the Bloch functions, which are periodic 
in the unit cell. Inserting these into Eqs. (O, one obtains 

^KO - ^r5]5]^|z..,.,(o)p|^,,;,,(o)p 

X Tr [Gi(fei, w„)ir2Gj(fc2, ~{ujn + Vm)) 

xfij(fci,fc2;Wn,^'m)] , (5) 



N{u:) = ^^„,.Re 

i 

M{uj) = ^7V„,,Re 



Here 



(8) 



(9) 



where Gi{k,ujn) are the impurity-averaged Green's func- 
tions given by 



Gi(fc,a;„) 



^^nZi(tJ„)To -I- ^i,fcT3 -I- 4)i{uJn)Tl 

c.2zf(^„) + eL + '/'?K) 



(6) 



Here Zi(ijJn) and (j)i{uin) are the renormalization function 
and the pairing self-energy, respectively, for the ith band. 

In Ref. we argued that the contribution to the 
vertex functions F,y(fei, fc2; i^m) from the electron- 
phonon interaction can be ignored while the effect of the 
Coulomb interaction drops out from the ratio Rg/ Rn- In 
calculating the contibution to F.^ from the impurity scat- 
tering in the conserving approximation p3 | one considers 



(10) 



are the local densities of states at the nuclear site in the 
normal state, with Np^i the Fermi level density of states 
in ith band, iV„ = Nn,i+Nn^2, and Aj(cj) = 0j(Li;)/Zi(cj) 
is the gap function in ith band. 

The gap functions are obtained by solving the finite 
temperature Eliashberg equations on the real axis |3Cl|| . 
which include the electron-phonon interaction, screened 
Coulomb repulsion and both the interband and intra- 
band impurity scattering described by the self-consistent 
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second Born approximation: 



(11) 



Ul' 

xK+^,j{uj, ~uj') - fi*j{ujc) tanh — 

"fzj 1 



(12) 
(13) 



+KT%{cj,cj')+K^%{cj,~u;')] , (14) 



where 



+ 00 



1 



± 



1 



uj' + LU + Q + iO+ uj' -LU + n-iO+ 



(15) 



1 



± 



1 



uj' + l^J + Q. + iO+ uj' - UJ + Q. - 



(16) 



Here Q;^i^y(r2) and fi*j{ujc) are intraband and inter- 
band electron-phonon coupling functions and Coulomb 
repulsion parameters for the cutoff lUc, respectively. 
The impurity scattering rates are defined by 7y/2 



■KNF,j{Q)\Vij\'^ where 



is the concentration of 
impurities and Vij is the Fermi surface averaged matrix 
element of the change in the lattice potential caused by 
an impurity between the states in the bands i and j. 

It is easy to see that the intraband scattering rates 
7m drop out from the equations for the gap functions 
Ai(aj) = (f)i{Lo) / Ziiu) which are obtained from Eqs. Ijllt 
I14|l . and only the interband impurity scattering affects 
the gap functions and hence Rs/Rn, Eqs- H7I9|) . We 
should point out that the Eqs. pill4|l are the same as 
the strong coupling equations for the McMillan tunnel- 
ing model of the superconducting proximity effect [sij 
and were first solved numerically at zero temperature by 
Zarate and Carbotte [s^l over twenty years ago. 
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FIG. 1: The ratio Rs/Rn in the case when the relaxation is 
dominated by the lower-gap band (A''„,ct=0) for several values 
of the interband scattering rate 7^0- = 7 in the units of tran- 
sition temperature of the clean system Tco- In (a) Rs/Rn is 
plotted as a function of temperature in units of Tco- Note that 
for 7/T'cO=10 the results do not extend all the way to the tran- 
sition temperature Tc because of the poor convergence of the 
real-axis Eliashberg equations near Tc for such a high value 
of 7. In (b) Rs/Rn is plotted as a function of T /Tc, where Tc 
is the transition temperature of disordered system, in order 
to illustrate the change in the relative width and position of 
the coherence peaks with increasing 7. 



III. RESULTS 

In order to examine the effect of the interband impu- 
rity scattering on the NMR relaxation rate of a singlet 
two-band superconductor we have calculated Rs/Rn for 
the interaction parameters of MgB2 . The four electron- 
phonon coupling functions cP'Fij (fi), i, j = cr, tt, were cal- 
culated by Golubov et al. 5] and the Coulomb repulsion 
parameters fj,*j{ujc) were determined in p^l27j| based on 
the screened Coulomb repulsion parameters of MgB2 cal- 
culated by Choi et al. 33]. Since j^^/j^^^ — Np^j^/Np^cr 
there is only one independent interband scattering rate 
parameter and we choose 7^0- = 7. Our representative 
results are shown in Figs.n]and[21 Figurenis our theoret- 
ical prediction for the relaxation rates on ^^Mg nuclei in 
MgB2, for which the dominant relaxation mechanism is 
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FIG. 2: The ratio Rs/Rn for iV„,^/iV„,^=0.45 for several val- 
ues of the interband scattering rate 7,ro- = 7 in the units of 
the transition temperature of the clean system Tco- In (a) 
Rs/Rn is plotted as a function of temperature in units of Tco- 
In (b) Rs/Rn is plotted as a function of T /Tc, where Tc is the 
transition temperature of disordered system. 



the Fermi contact interaction |3J,|35j and the electronic 
structure calculations [s^ give A^n.cr ~0 for the Mg-site. 
To the best of our knowledge there are no measurements 
of l/Ti in the superconducting state on ^^Mg in MgB2, 
presumably because of its small magnetic moment and a 
low natural abundance, but the experiments performed 
in [33, indicate that it is possible to measure ^^i? be- 
low the superconducting transition temperature. Such 
measurements would be highly desirable since our the- 
ory is quantitatively correct for ^^Mg nucleus. The broad 
peaks in Fig. [H between O.STc and 0.6Tc for 0< 7 <2Tc 
are analogous to the broad peaks found in the microwave 
conductivity of MgB2 js^ in the same temperature range. 
In the hest fits to the data were obtained by assum- 
ing that the 7r-band gives the dominant contribution to 
the microwave conductivity and our results in Fig.^Eilso 
give only the 7r-band contribution to Rs/Rn {Nn,a = 0). 
Since both NMR relaxation rate and the microwave con- 
ductivity have the same coherence factors in the single 
band case |40j the similarity between our prediction for 
^^R and the results obtained in [s^ is not accidental. 

In Fig. 121 we present our results for Nn-nlNn,a = 0.45, 
which would correspond to "'^^B nucleus j3(|| in MgB2 if 
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FIG. 3: The densities of states, Re [u/ \J lo^ — Ao-,7r(cL))], in 
the a and tt bands at (a) T = and (b) T = O.gSTc, for the 
impurity interband scattering rates 7 = O.lTco and 7 = Tco- 



the dominant relaxation mechanism were the Fermi con- 
tact interaction. However, the local-density approxima- 
tion (LDA) calculations |3J, |33 have found that at the 
^^B nucleus the most significant contribution to the re- 
laxation comes from the interaction with the electronic 
orbital part of the hyperfine field. Hence, our results in 
Fig. 121 should not be compared directly with the experi- 
mental results for ^^R. While the predictions in |33.l35l | 
were confirmed by experiments |37l IsM l4lj| in the nor- 
mal state, until recently |42l | nothing was known theo- 
retically about the temperature dependence of 1 /Ti in a 
superconductor in which the orbital part of the hyper- 
fine field dominates the NMR relaxation. Our prelimi- 
nary results j42i] for a single band superconductor indi- 
cate that in such a case the temperature dependence of 
Rs/Rn is given by the standard expressions obtained for 
the Fermi contact interaction 30, 40J provided that 7 is 
much greater than Tc- In our treatment ^| of the or- 
bital contribution to Rs/Rn the extended nature of the 
single electron states participating in the formation of 
Cooper pairs pla yed the key role, and it is not clear how 
the treatment of |34Ll3q | in which a few localized orbitals 
are responsible for the orbital part of the hyperfine field 
could be extended to the superconducting state. 

It is apparent from Figs. and [21 that as the inter- 
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band impurity scattering rate increases initially from 
to about 2Tco the coherence peak in the NMR relaxation 
rate is reduced and moved closer to the Tc- In the case 
of ^^i? the peak is also broadened (see Fig.QjD) while the 
shoulder in Rs/Rn for Nn,^/Nn,a = 0.45 in Fig. for 
7=0 resulting from the 7r-band contribution is rapidly 
reduced, Fig. ^ As we have anticipated in Sec. ^ these 
changes in Rs/Rn at low values of 7 are a direct con- 
sequence of the changes in both normal and anomalous 
densities of states in the superconducting state with in- 
creasing impurity scattering, Eqs. (|7I9|I . In Fig. |3| we 
show the partial densities of states in the two bands at 
T = and at T just below the T^, for two lowest values 
of the interband impurity scattering parameter 7 from 
Figs, n and El Clearly, in the low 7 regime increasing 
interband impurity scattering leads to a reduction and 
smearing of the partial densities of states both at low and 
high temperatures. As a result the the Hebel-Slichter co- 
herence peaks are reduced in size and in the case of a 
single-band contribution. Fig. ^ the peak is also broad- 
ened. 
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nearly equal, resulting in nearly equal normal and anoma- 
lous densities of states in both bands. As a result 
one expects the system to display l/Ti characteristic 
of a single-band superconductor with pronounced Hebel- 
Slichter peak, barring extremely strong elcctron-phonon 
coupling H^ll^. Indeed, as illustrated in Fig. ^ for 7 
= lOTco, the densities of states in the two bands are 
very similar, in particular at high temperature where 
the damping term in Eq. associated with the inter- 
band impurity scattering is reduced by the smaller gaps 
at high temperatures. The corresponding Rs/Rn shown 
in Figs, n^a) and |2Ia) have well pronounced coherence 
peaks similar in shape to what one would expect for a 
single-band s-wave superconductor. 



IV. CONCLUSIONS 

We calculated the NMR relaxation rate 1/Ti in an s- 
wave two-band superconductor with impurities, assum- 
ing that the relaxation of nuclear spins is controlled by 
the Fermi contact interaction with the band electrons. 
Using the interaction parameters of MgB2 0, we 
found that for low interband impurity scattering rates 
disorder suppresses Hebel-Slichter coherence peak as a 
result of the smearing of the densities of states in the 
two bands with two different gaps. For high scatter- 
ing rates, as the gap functions in the two bands become 
nearly equal, the system behaves as a single-band super- 
conductor with a well-developed coherence peak below 
Tc, as appropriate for a medium electron-phonon cou- 
pling parameter A ~1 and a small ratio Tc/^max, where 
^max is the maximum phonon frequency. However, one 
should keep in mind that in the limit of strong disorder, 
such as found in the samples of MgB2 irradiated with 
high neutron fluences JJj _19, _20, 21J, our treatment of 
disorder is not sufficient to quantitatively describe all ex- 
perimental results. Indeed recently observed [3 l2(l l2ll | 
reduction of the smaller 7r-band gap with increasing dis- 
order in neutron irradiated samples of MgB2 can never 
be reproduced within the self-consistent second Born ap- 
proximation treatment of impurity scattering used here 
and elsewhere in the literature on two-band superconduc- 
tivity. For a highly disordered system one would likely 
have to consider changes in the electron-phonon interac- 
tion and the normal state densities of states to account 
for the reduction of the gap in three-dimensional 7r-band. 



FIG. 4: The densities of states. Re [ui/y^uj'^ — Ao-,7r(i^)], in the 
a and n bands at (a) T = 2 K and (b) T = 28 K for the highest 
impurity interband scattering rates 7 = lOTco considered in 
this work. 

At high vales of 7 (the Anderson limit) one expects 
that the gap functions in the two bands will become 
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